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摘要 ”Fe-Cr 合金 作为 聚变 反应 堆 候选 材料 中 低 活 化 铁 素 体 钢 ( 马 氏 体 钢 ) 的 主体 成 分 ， 在 实际 工 况 下 ， 聚 
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变 反 应 产生 的 氨 会 进入 材料 中 ， 产 生 微观 缺陷 ， 从 而 引起 辐 照 损伤 ， 影 响 材料 的 抗 辐 照 能 力 及 其 他 物理 化 


学 性 能 。 本 文通 过 第 一 性 原理 计算 ， 模 拟 得 出 反应 堆 主 体 材料 Fe-Cr 合金 的 界面 ， 并 在 不 同 的 间隙 位 置 和 


BUBBA He 原子 ， 结 构 优 化 后 计算 形成 能 ， 得 出 He 原子 在 不 同位 置 的 稳定 性 。 与 此 同时 ， 通 过 对 体 


积 的 计算 得 出 He 在 靠近 界面 的 Fe 八 面体 间隙 位 时 会 引起 较 大 的 晶 格 畸变 ， 这 是 He、Fe、Cr 三 种 原子 在 
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特定 位 置 相互 杂 化 影响 的 


关键 词 Fe-Cr 合金 ”第 一 性 原理 Beek FM 


ur 
m 


First principles study on helium doping at the interface of Fe Cr alloy 
LIAN Xiang-Yu!? 
1(Anhui Nuclear Exploration Technology Central Institute, Wuhu, 241000, Anhui) 
2(Radioactive Resources and Environmental Engineering Research Center of Anhui Provincial Bureau of Geology and Mineral 


Resources, Wuhu, 241000, Anhui) 


Abstract: As the main component of low-activation ferrite/martensite steel in candidate materials for fusion 
reactor, Fe-Cr alloy will enter the material and produce micro-defects under actual working conditions, which will 
cause radiation damage and affect the radiation resistance and other physical and chemical properties of the 
material. In this paper, the interface of Fe-Cr alloy, which is the main material of reactor, is simulated by 
first-principles calculation. He atoms are doped in different gap positions and substitution positions. The formation 
energy is calculated after structural optimization, and the stability of He atoms in different positions is obtained. At 
the same time, the calculation of volume shows that the lattice aberration of He at the interstitial position of Fe 


octahedron near the interface is caused by the hybridization of He, Fe and Cr atoms at specific positions. 
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随 着 社会 的 发 展 和 人 类 的 不 断 进步 , 能 源 短缺 问题 正 逐 步 受到 重视 , 核能 因此 成 为 了 一 
种 新 的 重要 能 源 。 随 着 反应 堆 的 建立 和 及 更 新 换代 ， 核 反应 堆 材 料 的 研究 也 变 得 日 益 
Fe-Cr 合金 作为 聚变 反应 堆 候 选材 料 中 低 活化 铁 素 体 / 马 氏 体 钢 的 主体 成 分 ， 长 期 以 来 受到 
广大 学 者 的 关注 与 研究 "。 国 内 外 低 活 化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 如 表 1 所 示 。 在 实际 
工作 条 件 下 , 材料 在 高 能 中 子 的 辐 照 下 会 产生 大 量 的 微观 缺陷 , 其 中 包括 间隙 型 ， 空 位 型 ”， 
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位 错 号 等 ， 微 观 缺 陷 及 其 引起 的 辐 照 损伤 会 对 材料 的 物理 化 学 性 能 产生 重要 的 影响 。 因 此 ， 
本 文选 择 Fe-Cr 合金 为 研究 材料 。 
表 1 国内 外 低 活 化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 [4-9] 

Tab. 1 Main components of low activated ferrite / martensitic steel all around the world'*! 
国家 型 号 主要 合金 元 素 微量 合金 元 素 
欧洲 EUROFER97 Fe (288) ,9Cr W, V, Ta 
美国 9Cr2WVTa Fe 基体 ) ,9Cr 2W, V, Ta 
日 本 F82H Fe (基体) ,7.5Cr 2W, V, Ta 
日 美 联 合 JLF-1 Fe 基体) ,9Cr 2W, V, Ta 
日 美 联 合 JLF-2 Fe 基体) ,9Cr 2W, V, Ta 
日 美 联 合 JLF-3 Fe (基体 ) ,7Cr 2W，V，Ta 
中 国 CLAM Fe (基体 〉 ,9Cr 1.5W, V, Ta 
中 国 CLF-1 Fe (基体 ) ,8Cr 1.5W, M, V, Ta 
中 国 SCRAM Fe《〈 基 体 ) ,9Cr 2W, V, Ni, Ti 


近 几 年 来 ， 基 于 Fe-Cr 合金 中 的 空位 ， 位 错 以 及 Cr 析出 相 ， 大 量 的 不 同 尺度 的 研究 
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直 在 进行 中 ， 有 文献 表明 ， 相 邻 位 置 的 Cr 原子 相互 吸引 ,， T 
相互 作用 是 较 弱 的 ”。 与 此 同时 ， 与 别 的 相关 元 素 ,，( 如 M 
有 效 的 禁止 空洞 和 位 错 环 的 扩张 ”在 空冷 济 火 的 Fe-Cr 样品 中 , 当 退 火 温度 上 升 至 773 K 
时 ， 少 量 的 位 错开 始 消失 ， 同 时 ，Cr 原子 在 900 K 时 发 9 


而 Fe-Cr 合金 


Cr 原子 与 空位 的 


n) 相 比 ，Fe 中 添加 的 Cr 原子 会 


ERSE, TE 1073 时 析出 物 溶解 "。 


在 这 些 研究 基础 上 , 通过 离子 注入 以 及 电化 学 充 入 等 方式 , 同时 展开 了 大 量 的 Fe-Cr 模型 钢 


六 此 
中 和 毛 行 为 的 研究 。 实 际 应 用 中 ， 核 反应 产生 的 He 原子 容易 进入 材料 


聚集 在 氨 泡 上 ， 即 氨 泡 的 长 大 。 而 对 于 材料 内 部 的 气泡 ， 


而 形成 非 本 征 缺 陷 。 


两 个 及 以 上 的 He 原子 与 气泡 结合 形成 氨 泡 ， 接 下 来 周围 的 空位 以 及 很 容易 迁移 的 氨 原 子 会 
' 情 况 是 辐 照 下 残留 的 气体 与 贫 
原子 区 和 过 饱和 空位 浓度 之 间 的 相互 作用 , 另 一 种 情况 是 妨 变 产物 与 气泡 相 结合 。 并 且 在 辐 


照 下 ， 材 料 中 很 容易 产生 Frenkel 缺陷 对 。 如 果 气 泡 增 大 ， 即 较 多 的 空位 进入 气泡 ， 就 意味 


量 He 的 扩散 是 氨 泡 成 核 和 生长 所 必须 的 ， 单 个 He 原子 


着 与 之 对 应 的 间隙 位 He 原子 进入 位 错 和 位 错 坏 ， 从 而 引起 品格 畏 变 ， 导 致 材料 的 肿胀 。 大 
易于 迁移 > 


且 与 空位 和 其 他 He 原子 


聚集 形成 He-V《〈 空 位 ) Bi, TRI He- 空 位 团 簇 与 其 他 的 He 原子 ，He 


7 团 簇 结合 ， 最 终 产 生 纳米 尺寸 的 氨 泡 。 


通过 基于 密度 泛 函 理 论 的 第 一 性 原理 计算 ， 可 以 高 
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移行 为 的 研究 提供 理论 基础 。 


精度 的 计生 
掉 不 同位 置 的 稳定 性 ， 从 理论 上 可 以 得 出 He 对 Fe-Cr 界面 唱 格 参数 的 影响 ， 为 界面 处 氨 迁 


kE px 


WR» 空位 以 及 He- 


得 出 He 原子 在 Fe-Cr 界 


1. 计算 方法 
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第 一 性 原理 是 


AHH 


任何 经 验 参 数 , 仅 根 据 粒 子 与 粒子 间 的 相互 作用 原理 ， 
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的 收敛 条 件 够 高 , 理论 上 能 精确 
可 分 为 基于 密度 泛 函 理论 的 计算 和 从 头 ? 


算 为 基础 的 。 


文件 包括 : 


INCAR, POSCAR, 


K 点 的 选择 ， 即 网 格 的 选择 。POTCAR 是 材料 计生 


计算 材料 各 方面 


目前 ， 第 一 性 原理 计算 中 ， 最 常用 


的 计算 


8 材料 的 各 


POTCAR 和 KPOI 


料 结构 的 特点 ， 选 用 不 同 的 帮 势 进行 计算 。 
坐标 。INCAR 则 是 VASP 软件 中 的 控制 文人 
的 收敛 条 件 ， 选 取 的 模型 方法 ， 以 及 计算 后 的 输入 文件 及 内 容 等 。 


Sits 


2. 结构 优化 


首先 , 间 际 位 置 包 括 四 面体 间 际 位 和 八 面 
结构 都 属于 在 体 心 立方 结构 (bcc)。bcc 唱 体 


TS。 其 中 KPOINTS 决定 了 计算 过 程 
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通过 量 


性 质 的 方法 。 也 就 是 说 ， 只 要 计算 资源 足够 ， 设 置 
' 物 理化 学 相关 性 质 。 广义 的 第 一 性 


原理 


十 算 ， 其 中 从 头 计算 是 以 Hartree-Fock 的 自治 场 计 
的 软件 就 是 维也纳 从 头 计算 模拟 包 VASP) ,输入 


布 里 渊 区 


所 选用 的 厦 势 ,不 同 的 材料 需要 根据 其 
POSCAR 决定 了 原子 的 品格 参数 ， 体 积 


材 
小 及 其 


EF， 也 就 是 所 谓 的 程序 执行 者 ， 不 同 参 数 的 设置 决 


体 间 隙 位 。 无 论 是 单 品 Fe 还 是 站 


HHCY > 


其 原子 


间隙 位 ， 因 此 一 般 情 况 下 ,间隙 原子 应 首先 


5 据 


中 计算 模型 为 2X2X4 的 Fe-Cr 超 唱 胞 ， 计 和 


掉 体 间隙 位 具有 不 对 称 性 , 如果 八 面体 间隙 位 被 间隙 原子 占据 ， 只 会 引起 吕 
近 (a/2) 的 两 个 原子 明显 的 偏离 平衡 位 置 ， 
显 偏 移 ， 因此 不 同体 心 立方 结构 金属 
性 原理 ， 通 过 VASP 软 件 ， 对 He 原子 和 空位 在 Fe-Cr 合 金 界 面 不 同位 置 处 的 结构 进行 优化 ， 其 


的 间 际 原子 占据 的 


， 四 面体 间隙 位 的 空间 体积 大 小 大 于 八 面 体 
四 面体 间隙 位 。 然 而 ， 体 心 立方 结构 中 的 八 


E 离 间隙 中 心 位 最 


间隙 原子 周围 剩 下 的 四 个 原子 则 并 不 会 发 生 明 


过 程 采 用 交换 关联 泛 函 至 


Hid (GGA) 和 


又 
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位 置 是 不 一 样 的 。 本 论文 基于 第 一 


加 平面 


波 方法 (PAW)“"#。 根 据 实验 前 的 调试 ，K 点 选择 Monhkorst-Pack 方 法 ， 结 构 优 化 时 设置 为 
6X6X4。 为 了 确保 品格 参数 及 其 形状 大 小 更 符合 实际 情况 ， 在 计算 之 前 ， 蝇 胞 形状 和 各 原 
子 的 位 置 得 到 完全 的 优化 ， 其 中 ISIF 设 置 为 3， 平 面 波 截断 能 (ENCUT) 设置 为 400 ev， 收敛 


条 件 设 置 为 每 个 原子 所 受到 的 力 小 于 0. 001 eV, 优化 后 得 到 合 弄 


并 在 此 基础 上 进行 静 


人 确 


态 计算 ， 


的 原子 分 布 模型 ， 如 图 1 所 示 。 


计算 过 程 


的 原子 排 布 结构 和 品格 参数 ， 


之 前 的 部 分 参数 不 变 ， 设 置 IBRION=-1， 得 出 较 精 


图 1 Fe-Cr 合金 界面 的 原子 模型 (1,2,3,4 分 别 代表 平行 于 xoy 平面 的 切面 ， 黄 色 代表 Fe 原子 ， 蓝 色 代 表 Cr 


原子 ) 
Fig. 1 Atomic model of Fe-Cr alloy interface (1, 2, 3, 4 represent the sections parallel to the xoy plane, respectively, 


yellow represents the Fe atom and blue represents the Cr atom. ) 


RIF, HERK Cr 原子 ， 黄 色 代表 Fe 原子 ,切线 1、2、3 和 4 分 别 表示 He 原子 所 


在 的 平面 。 另 外 ，T (Tetrahedron) 表示 He 原子 位 于 四 面体 间隙 位 ，0 (0ctahedron) 
表示 He 原子 位 于 八 面体 间隙 位 、S CSupersede) 表示 He 原子 位 于 蔡 代 位 ， 为 了 进行 对 
比 ， 另 加 入 空位 : V Vacancy) 表示 空位 在 替代 位 《图 中 某 个 Fe 或 Cr 原子 消失 ， 其 原 


所 在 的 位 置 没 有 任何 原子 )。 比 如 : T3 就 表示 He 原子 位 于 3 号 切面 上 ， 并 处 于 由 一 个 Cr 
原子 和 三 个 Fe 原子 所 构成 的 四 面体 的 中 心 。 在 结构 优化 的 基础 上 ， 分 别 把 He 原子 添加 
到 以 上 位 置 处 ， 再 对 其 进行 结构 优化 ， 优 化 参数 与 之 前 相同 ， 结 构 优 化 完成 后 进行 静态 
计算 。 为 了 得 出 He 原子 在 各 个 不 同位 置 处 的 稳定 性 ， 并 利用 以 下 公式 计算 其 形成 能 ”: 


FFormation = Epe, cr. He 一 Ere cr 7 Ene (1) 


EFormation = Ere. Cr. He 一 Ere. Cr EHe + Esingle Fe/Cr (2) 


HAA, Er c te KRÉERA He 原子 的 Fe-Cr 界面 体系 的 总 能 量 ，B. 表示 Fe-Cr 界 


面体 系 的 总 能 量 ，Bi 表示 单个 He 原子 在 一 个 很 大 空间 内 的 总 能 量 ”。E,iworwe 表示 被 营 


代 的 单个 Fe 或 Cr 原子 的 静止 能 量 。 在 经 过 相关 计算 过 后 ,得 出 来 的 He 在 Fe-Cr 界面 各 


个 位 置 的 形成 能 如 图 2 所 示 。 
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FeCrV1 V2 V3 V4 S1 S2 S3 S4 01 02 O3 O4 T1 T2 T3 T4 - 
He 和 空位 所 在 的 不 同位 置 
图 2 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 形成 能 


Fig. 2 Formation Energy of He and Vacancies at Diffe rent Positions of Fe-Cr Alloy Interface 


3. 形成 能 


MEI VI, V2, V3 及 V4 的 形成 能 可 以 得 出 ， 当 Fe-Cr 界面 中 少 了 一 个 原子 ， 导 致 
单 空位 产生 时 ， 其 结构 是 不 稳定 的 。 相 比 于 Cr 原子 处 出 现 空位 ， 空 位 更 倾向 于 在 远离 界 
Hi Fe 原子 位 置 形 成 (V4 位 置 表示 远离 界面 的 Fe 原子 位 置 ， 其 形成 能 相对 最 小 ， 稳 定 
性 也 相对 最 高 )。 换 句 话 来 讲 ， 当 Fe-Cr 合金 中 出 现 单 空位 时 ， 单 空位 更 倾向 于 在 远离 界 
ANY Fe 原子 位 置 处 形成 ， 即 四 周 都 是 Fe 原子 的 位 置 。 这 为 接 下 来 研究 Fe 基 合 金 核反应 
堆 材料 中 空位 的 形成 位 置 ， 提 供 了 理论 依据 。 


根据 图 中 He 原子 在 不 同位 置 (S1，S$2，S3，S4，01，02，03，04，T1，T2，T3，T4) 


的 形成 能 可 以 得 出 ，He 在 四 面体 间隙 位 最 稳定 ， 其 次 是 He 在 八 面体 间隙 位 ，He 在 替代 
位 相对 最 不 稳定 。 在 这 之 中 ，Tl 位 置 的 形成 能 最 低 ，S1 位 置 的 形成 能 最 高 ， 可 以 得 出 ， 
当 一 个 He 原子 进入 Fe-Cr 界面 处 时 ， 会 首先 倾向 位 于 远离 界面 的 Cr 四 面体 间隙 位 中 ， 
多 个 He 原子 进入 后 会 对 其 他 的 四 面体 间隙 位 进行 填充 ， 直 至 靠近 界面 的 Fe 四 面体 间 际 
位 都 填 满 后 ， 才 会 进入 八 面体 间隙 位 中 ， 一 直到 八 面体 间隙 位 也 填 满 ， 过 多 的 He 原子 就 
会 开始 蔡 代 原 有 的 Cr 原子 与 Fe 原子 ， 取 代 他 们 的 位 置 ， 形 成 氨 泡 或 所 空位 复合 体 ， 对 


核反应 堆 材 料 的 抗 辐 照 能 力 等 产生 较 大 影响 。 而 被 取代 的 Cr 原子 与 Fe 原子 ， 由 于 脱离 
了 原 有 的 品格 位 置 ， 会 因为 较 大 的 原子 半径 ， 对 原 有 晶 格 参数 产生 较 大 影响 ， 在 材料 中 
形成 微观 缺陷 ， 如 空位 、 空 位 团 、 位 错 、 位 错 环 等 ， 进 一 步 影 响 材 料 的 物理 化 学 性 质 。 


4, 体积 变化 


为 了 更 直观 的 得 出 He 在 Fe-Cr 界面 不 同位 置 的 影响 ， 第 一 性 原理 计算 得 出 He 和 空 
位 在 Fe-Cr 界面 不 同位 置 处 的 体积 如 图 3 所 示 。 
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图 3 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 体积 


Fig. 3 Volume of He and Vacancies at Different Positions of Fe-Cr Alloy Interface 


显而易见 ， 相 比 于 没有 任何 缺陷 的 Fe-Cr 合金 界面 , 无 论 是 He 在 替代 位 还 是 空位 在 
蔡 代 位 ， 都 会 使 得 Fe-Cr 界面 的 体积 相对 减 小 ， 且 单 空位 的 出 现 ， 更 会 明显 降低 Fe-Cr 
界面 的 体积 。 当 He 原子 出 现在 间 隐 位 时 ， 不 论 是 四 面体 间隙 位 还 是 八 面体 间隙 位 ， 都 会 
明显 的 增加 Fe-Cr 界面 的 体积 。 值 得 注意 的 是 ， 当 He 原子 出 现在 Fe-Cr 界面 的 不 同位 置 
时 CO1, 02, 03, 04, T1, T2, T3, T4), T3 位 置 的 体积 比 其 他 位 置 的 体积 都 大 ， 且 体 
职 的 变化 近乎 是 别 的 位 置 的 两 倍 。 计算 结果 表明 ， 当 He 原子 位 于 靠近 界面 的 Fe 八 
自体 间隙 位 时 ， 会 引起 整个 界面 出 现 较 大 的 晶 格 畸变 ， 这 种 畸变 会 一 定 程度 上 改变 材料 
的 形状 ， 从 而 影响 材料 的 物理 化 学 性 质 。 
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针对 于 这 一 位 置 的 特殊 现象 , 我 们 继续 对 He 在 此 位 置 的 结构 模型 进行 了 多 次 精度 更 
高 的 优化 ， 其 优化 结果 都 在 358. 94—359. 2 让 之 间 (1A=10"m) 。 与 此 同时 ， 使 用 p4vasp 
软件 得 出 了 它 的 态 密度 图 以 及 包含 不 同 电子 轨道 的 分 波 态 密度 图 ， 并 与 其 他 八 面体 间 辽 
fi COL, 02,04) 的 电子 态 密度 图 进行 了 比较 。 
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Al 4 He 在 Fe-Cr 合金 界面 不 同 八 面体 间隙 位 01，02，03，O04 的 电子 态 密度 图 ， 横 坐标 表示 能 量 区 间 ， 纵 


| 之 


标 表示 电子 态 密度 


(03 中 红 绿 蓝 三 条 曲线 分 别 表示 Fe、Cr、He 原子 的 分 波 电子 态 密度 ) 


Fig. 4 Electronic Density of State Charts of He at Different Octahedral Gap Points O1, O2, 03, O4 at Fe-Cr Alloy 
Interface (The three curves of red, green and blue in O3 represent the fractional density of states of Fe, Cr and He 


atoms, respectively ) 


电子 态 密度 表示 单位 能 量 范围 内 所 允许 的 电子 数 ， 也 就 是 说 电子 在 茶 一 能 量 范围 的 


分 布 情况 。 因 为 原子 轨道 主要 是 以 能 量 的 高 低 去 划分 的 ， 所 以 电子 态 密度 图 能 反应 出 电 
子 在 各 个 轨道 的 分 布 情况 ， 反 映 出 原子 与 原子 之 间 的 相互 作用 情况 。 

从 图 中 可 以 看 出 ， 其 他 八 面体 间隙 位 (01，02，04) 在 能 量 -20eV 一 -18eV 处 有 一 
个 波峰 ， 而 He 在 靠近 界面 的 Fe 八 面体 间隙 位 (03) 时 的 波峰 出 
之 间 。 波 峰 的 偏 移 表 示 He 原子 在 此 位 置 会 与 周围 的 Cr 原子 和 Fe 原子 有 特殊 的 相互 作用 ， 
该 峰 位 不 只 是 一 个 原子 作用 的 结果 ， 而 是 Fe、Cr、He 三 种 原子 所 代表 的 红 绿 蓝 三 条 线 都 
有 波峰 产生 ， 说 明 He 原子 的 电子 轨道 与 Fe 原子 和 Cr 原子 的 电子 轨道 出 现 了 杂 化 现象 。 
该 峰 是 由 Fe, Cr, He 所 共同 组 成 的 杂 化 峰 ， 是 Fe. Cr. He 三 原子 之 间 一 种 特殊 的 相互 
作用 。 但 是 相 比 于 其 他 八 面体 间隙 位 来 说 不 稳定 ，He 原子 很 容易 迁移 到 别 的 八 面体 间隙 
位 以 及 四 面体 间隙 位 。 


现在 能 量 -18eV 一 -16eV 


第 一 性 原理 计算 得 出 ， 当 He 原子 进入 Fe-Cr 合金 界面 时 ， 首 先 倾向 于 向 远离 界面 位 
置 的 Cr 四 面体 间隙 位 置 处 移动 , 当 这 一 位 置 被 完全 占据 , 之 后 依次 是 其 他 四 面体 间隙 位 ， 
八 面 体 间 阶 位 ， 蔡 代位 。 这 一 结论 有 助 于 进一步 研究 核反应 堆 材 料 中 的 氮 迁 移行 为 ， 从 
而 研究 材料 的 抗 辐 照 能 力 。 并 且 ， 当 He 原子 进入 靠近 Fe-Cr 合金 界面 的 Fe 八 面体 间 际 
位 时 ， 会 引起 较 大 的 晶 格 畸变 ， 这 是 Hey Fey Cr 三 种 原子 在 特定 位 置 相互 杂 化 影响 的 结 
果 。 这 使 得 界面 的 原子 排列 发 生 较 大 改变 ， 从 而 影响 材料 的 微观 缺陷 结构 。 
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